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Abstract 

The mechanism of ferromagnetism stability in sandwich clusters V„(C6He) n+ i has been studied 
by first-principles calculation and model analysis. It is found that each of the three types of bonds 
between V and benzene (Bz) plays different roles. V 3d 2 2 orbital, extending along the molecular 
axis, is weakly hybridized with Bz's HOMO-1 orbital to form the a-bond. It is quite localized and 
singly occupied, which contributes \\ib to the magnetic moment but little to the magnetic coupling 
of neighboring V magnetic moments. The in-plane d x 2_ y 2, d xy orbitals are hybridized with the 
LUMO of Bz and constitute the 5-bond. This hybridization is medium and crucial to the magnetic 
coupling though the 5 states have no net contribution to the total magnetic moment. d xz , d yz and 
HOMO of Bz form a quite strong 7r-bond to hold the molecular structure but they are inactive in 
magnetism because their energy levels are far away from the Fermi level. Based on the results of 
first-principles calculation, we point out that the ferromagnetism stability is closely related with 
the mechanism proposed by Kanamori and Terakura [J. Kanamori and K. Terakura, J. Phys. Soc. 
Jpn. 70, 1433 (2001)]. However, the presence of edge Bz's in the cluster introduces an important 
modification. A simple model is constructed to explain the essence of the physical picture. 
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I. INTRODUCTION 



Since the discovery of ferrocene Fe(C5H 5 ) 2 and determination of its structure, intensive re- 
search of the sandwich complexes where a metal atom or ion is sandwiched with two organic 
molecules (mostly aromatic molecules) has been performed in the field of organometallic 
chemistry- The growing interests in these complexes are mainly due to their potential usages 
in catalysis, novel magnetic and optical materials and other various applications. Among 
these compounds, the 3d-transition-metal-organic-molecule complexes are particularly un- 
der current active study since molecular magnets are considered as potential candidates for 
future applications in high-density information storage. Therefore, since the discovery of 
ferromagnetism in vanadium-benzene V n Bz m complexes,- 1 ^ many experimental and theo- 
retical studies^&£ have been performed in order to shed light on the mechanism of their 
ferromagnetism stability. Now, it is generally accepted that for small size V n Bz m complexes 
(n<4), the energetically preferred structures are linear multiple-decker sandwich-like with 
m = n — 1, n and n + 1. The Stern-Gerlach deflection experiments show that the mag- 
netic moment in V n Bz n+1 increases linearly with n.— Ab initio calculations have confirmed 
this^ though the magnitude of the measured magnetic moment is smaller than the theo- 
retical value by about factor two, which is attributed to the spin relaxation effects during 
the measurement.-^ Also, for the ideal infinite (VBz) n=00 chain, first-principles calculations 
predicted that the ferromagnetic (FM) state is more stable than the anti-ferromagnetic 
(AFM) state.— A mechanism of double exchange interaction has been proposed^^ to 
explain the ferromagnetism in this complex based on the electronic structures obtained with 
the density functional theory calculations within the local density approximation (LDA) or 
generalized gradient approximation (GGA). 

However, it is well known that LDA or GGA does not reliably describe the localized state 
in transition metal compounds, especially for these low-dimensional complexes. Including 
on-site Coulomb interaction U usually can improve the description and lead to a more real- 
istic picture. In this paper, by performing GGA+U calculations, we find that the electronic 
structure is significantly modified. A mechanism, originally proposed by Kanamori and 
Terakura,— is found to be responsible for the stability of FM state both in infinite and finite 
V-Bz sandwich complexes. A similar mechanism was also proposed by Yabushita et al.->^ 
for the ferromagnetism stability in V2BZ3. However, it is also found that the presence of 
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edge Bz introduces an important modification for finite complexes. To show the essence of 
the physical picture, a simple model is constructed and the related parameters are fitted 
from the first-principles calculations of n = 2 cases. The numerical model analysis is semi- 
quant it at ively consistent with the first-principles calculations for n—2, 3 and 4 cases. The 
model calculation is also used to analyze the size dependence and configuration dependence 
of the energy cost for the magnetic moment reversal. 



II. METHODOLOGY 

All the first-principles calculations were performed using our software package OpenMX,— 
which is based on a linear combination of localized pseudo-atomic orbital (LCPAO) method. 
The PAOs are generated by a confinement potential scheme^ 4 - with a cutoff radius of 5.0 a.u. 
for hydrogen and carbon and 6.5 a.u. for vanadium, respectively. In the pseudopotential 
generation, the semicore 3s and 3p states of V atom were included as valence states. The 
exchange correlation energy functional within GGA as parameterized by Perdew, Burke and 
Ernzerhof (PBE)^ is used. Double-valence and polarization orbitals were included and the 
basis set completeness was checked by optimizing the structure of single Bz molecule and 
body centered cubic V metal and comparing with the experimental values or all electron full- 
potential calculations. The real-space grid technique^ was used with an energy cutoff of 250 
Ry in numerical integrations and in the solution of the Poisson equation using the fast Fourier 
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transformation algorithm. Hubbard U is included by the approach proposed in Ref. 
improve the description of localized states like V 3d orbitals. All the geometrical structures 
of V n Bz n+ i complexes are fully optimized under each individual magnetic configuration until 
the forces are less than 1.0 x 10 -4 a.u. 



III. RESULTS AND DISCUSSIONS 
A. Ideal infinite (VBz) n=oc chain 

To understand the mechanism of ferromagnetism stability in V-Bz multiple- decker sand- 
wich complex, we start by studying the bonding nature between Bz and V. In general, the 
molecular axis (z-axis) is defined as the line that passes through the metal and the center 
of gravity of Bz molecule and classify the valence orbitals in terms of their pseudo-angular 
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momenta around this axis. According to the symmetry, the six 7r-orbitals of benzene form 
one pa (HOMO-1), two prr (HOMO), two p5 (LUMO) and one pcf) (LUMO+1) orbitals. 
When sandwiched by two Bz molecules, the five 3d orbitals of the metal atom are classified 
to one da (d z a), two dn [d xz and d yz ), and two d5 (d xy and d x i_ y i) orbitals, and the 4s 
orbital is classified to an sa orbital.- 1 ^ To see the overall feature of the electronic structure, 
we performed first-principles calculations for both FM and AFM states of the ideal infinite 
V-Bz chain within GGA. The total and partial densities of states (DOS) for the FM state 
are plotted in Fig. 1, which are basically the same as other calculations.- 1 ^ 1 ^ By examining 
these DOSs and the symmetry of orbital wave-functions, one can easily find: 1) The coupling 
between pa and da is very weak, nearly non-interacting, making da localized and largely 
spin-split. 2) The prr and dn orbitals have strong hybridization, making both bonding and 
antibonding states far away from the Fermi-level, either totally occupied or totally empty. 
They are responsible to the stability of the molecular structure but not active in magnetic 
or electronic properties. 3) The p5-d5 hybridization is of medium strength among the three 
bonds. dS states are also spin split due to the Hund coupling with the spin polarized da 
states. V atom has 5 electrons, which occupy the da and doubly degenerate dS bands. 
Therefore, the number of holes existing in the overlapping da and d5 bands of the minority 
spin state is one. Due to this metallic feature, the double-exchange mechanism has been 
proposed to be the possible origin of ferromagnetism.- 1 ^ 1 ^ The FM state is more stable 
than the AFM state by 0.091 eV per VBz (see Table I). The total magnetic moment is just 
1 fiB per VBz unit with its integer value being consistent with the half- metallic feature of 
DOS. In fact, the local magnetic moment on V atoms is about 1.15 yUg. The total magnetic 
moment is reduced to 1 //g due to the negative magnetic moment on Bz being polarized 
through the p5-d5 hybridization.— 

However, for localized states like da and d5, LDA or GGA is not generally reliable. 
GGA+[7 could improve the description of the electronic structure. In the present calculation, 
U is a parameter but we have confirmed that the qualitative features are not affected by 
changing U from 2.0 eV to 4.0 eV. As a typical case, the results with [7=3.0 eV will be 
mostly shown below. The DOS of GGA+[7 calculation with [7=3.0 eV is shown in Fig.2a 
for the FM state. As expected, the spin-splitting in da and d8 is largely enhanced compared 
with that in GGA calculation. It is interesting to note that the broad band in the minority- 
spin state right below the Fermi level is mostly composed of the pS states, which are LUMO 
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FIG. 1: (Color online) DOS of a FM infinite V-Bz chain calculated within GGA. (a) Total DOS 
and the partial DOS of carbons in Bz molecules; (b) partial DOS for each of V 3d orbitals. 

in isolated Bz. Therefore,two dS electrons of V are formally transferred to p5 states of Bz. 
The spin splitting in dS states is now larger than the energy separation of d8 and p5 levels 
before the onset of the p — d hybridization. (More detailed discussion will be given below.) 
A similar calculation for the AFM state gives the DOS shown in Fig. 2b. Now, both the FM 
and AFM states are insulating and yet the stability of the FM state against the AFM state 
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FIG. 2: (Color online) DOS of an infinite V-Bz chain calculated within GGA+C/ (17=3.0 eV) for 
(a) FM and (b) AFM state. Upper panel of them is total and the partial DOS of carbons in Bz 
molecules, while the lower panel is partial DOS for each of V 3d orbitals. 
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FIG. 3: (Color online) A schematic picture describing the interactions between V db and intervening 
Bz p5 states in FM and AFM configurations. The thick vertical up arrow indicates the energy 
increasing direction and its left (right) side represents the majority- (minority-) spin channel. The 
horizontal solid (dashed) lines represent the energy levels of the corresponding state after (before) 
p — d hybridization. The hollow arrows indicates the direction and magnitude of movements for 
each state. Fermi level is labeled by the horizontal dash-dotted line. 
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is enhanced (see Table I). This implies that the mechanism of the ferromagnetism stability 
is not the double exchange. In the following, we further analyze the electronic structure in 
detail and elucidate the mechanism of the magnetic coupling between V magnetic moments. 

As da — pa hybridization is weak, da states does not make significant contribution to 
the magnetic coupling. Moreover, even if there is any contribution from them it will be 
AFM superexchange. The magnetic coupling is mostly governed by the d5 and pS states. 
Therefore, we will focus only on the S states. Figure 3 is a schematic diagram showing the 
basic feature of p — d hybridization in the 6 states. Due to the strong exchange splitting of 
d5 states, their majority-spin states are below the p5 states and fully occupied. On the other 
hand, the minority-spin d5 states are above the pS states and empty. In the FM configuration 
of the V magnetic moments, the majority-spin p5 states are strongly pushed up in energy by 
the p — d hybridization and vice versa in the minority-spin state. In the AFM configuration, 
on the other hand, the energy shift of p5 states is much reduced because the hybridization 
with the two neighboring dS states produces energy shifts with opposite sign. In addition 
to this, the energy shift is the same in both spin states. One additional comment is that 
the apparent exchange splitting of d5 states is enhanced by the p — d hybridization in both 
the FM and AFM states. Based on the observation of the energy shift scheme, we show a 
schematic DOS for S states in the FM state in Fig. 4. In the absence of p — d hybridization, 
the doubly degenerate p5 bands are half filled (Fig.4a). With the p — d hybridization which 
produces large energy shifts in the pS states (Fig. 3a), the DOS is modified to the one shown 
in Fig. 4b, which captures the essential feature of the S symmetry part of the DOS in Fig. 2a. 
Now, the electrons originally occupying the majority-spin p5 bands in Fig. 4a are transferred 
to the minority-spin pS bands as shown in Fig.4b, which produces a large energy gain in 
the FM state. A large negative spin polarization on Bz is also naturally explained. This is 
exactly the same story proposed by Kanamori and Terakura.— The importance of the spin 
polarization in the intervening Bz for the ferromagnetism stability was also pointed out by 
Yabushita et al. for V2BZ3 cluster.-^ 

As shown in Table I, the total energy difference between FM and AFM states increases 
as U value increases. Although the local magnetic moment on V ion increases with U, the 
negative polarization on Bz also increases and keeps the total magnetic moment per VBz 
unit the same as 1.0 fiB, which is produced by da states. 
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TABLE I: Magnetic moment on V (/is) and total energy difference between FM and AFM states 
(AE= Efm-^afm) calculated with different U values. 



U (eV) 


Infinite V-Bz Chain 


V 2 Bz 3 


Magnetic moment on V 


AE/V (eV) 


Magnetic moment on V 


AE/cluster (eV) 


0.0 


1.15 


-0.091 


1.10 


-0.003 


1.0 


1.35 


-0.113 


1.17 


-0.004 


2.0 


1.67 


-0.109 


1.29 


-0.004 


3.0 


2.03 


-0.184 


1.50 


-0.005 


4.0 


2.23 


-0.251 


1.82 


-0.007 




FIG. 4: (Color online) Schematic pictures describe the mechanism in which the p — d hybridization 
stabilize the FM states. The DOS (a) without and (b) with p — d hybridization is plotted. The 
energy gain due to the charge transfer is indicated. 

B. Finite V n Bz n +i clusters 

Now, we come to the study of finite size sandwich clusters V n Bz n+1 . For n=l, the spin- 
polarized VBz 2 with S=l/2 has lower energy than the non-magnetic one, which is consistent 
with the experimental measurement. For n—2, as shown in Table I the negative total energy 
difference between FM and AFM states indicates the FM state is more stable, being also 
consistent with experimental observations. However, the value of the energy difference is 
only several meV, much smaller than that in the infinite chain. Other calculations also give 
the energy difference less than 10 meV.-^ Before presenting further detailed analysis, the 
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following comment on Table I may be needed for the comparison of energies. The total 
energy difference is for one VBz unit in the infinite chain, while it is for V2BZ3 in the 
cluster. A justification of this seemingly strange comparison comes from an expectation 
that the magnetic coupling comes mostly from the energy change associated with the spin 
polarization of the intervening Bz between the V atoms. In this sense, both units have only 
one Bz shared by two neighboring V atoms. 

According to the physical picture described in the infinite chain, the magnetic coupling 
in V2BZ3 can be studied by a simple tight-binding model which takes account of the hy- 
bridization between p5 and dS states. For example the model Hamiltonian matrix for the 
AFM state is given as 

e p t' 

t' e d -A t 

t e p t • 

t e d + A t' 

t' e p 
where e p and e d are the on-site energies of p5 and d8 levels, respectively, and 2A is the 
exchange splitting. In order to see the role of the edge Bz explicitly, the hopping integral 
between d5 and p5 is expressed as t' for the edge Bz part and t for the middle Bz part. 

To remove the edge Bz's, we adopt a cyclic boundary condition in the (VBz) 2 cluster. 
The Hamiltonian, which is identical to the V point expression of an infinite chain, is written 
as 

Ep t t 

t e d -A t 

t E p t 

t t e d + A 

The eigenvalues of these Hamiltonians can be calculated easily. Assuming e p =e d =eo and 
taking the approximation that \t/A\ « 1.0 , we show the energy diagrams in Fig. 5. In 
the n = 2 case, there are 10 electrons introduced by two V atoms. Two of them occupy the 
da states and the other 8 electrons should occupy the dS and p8 states depicted here. Note 
that both dS and p5 states have double orbital degeneracy. Fermi level can be determined 
as shown in the energy diagrams. In order to understand the overall features of energetics 
among four different configurations shown in Fig. 5, we assume t' = t, which is reasonable 
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since the V-Bz distance is almost equally spaced. Then we note the following aspects in the 
energy diagram. 

1) As for the energy levels originating from d5 states, the sum of the occupied part of 
orbital energies is the same for all cases shown in Fig. 5. Note that the spin degeneracy has 
to be taken into account in the case of AFM state. Therefore, the total energy difference 
among the four cases shown in Fig. 5 comes entirely from the energy diagram associated with 
p5 states. 

2) In the case of FM state, even for the energy levels originating from pS states, the sum 
of occupied part of orbital energies is the same for both isolated cluster and cyclic boundary 
condition. 

3) In the case of AFM state with cyclic boundary condition, as there is no shift of energy 
levels of pS states, the FM state is definitely lower in energy by 8t 2 /A per cluster. 

4) In the case of AFM state for isolated cluster, the energy of one of the pS states is 
lowered due to the spin polarization at the edge Bz. Therefore the edge Bz contributes to 
the stabilization of AFM state of the isolated cluster. The FM state is still lower in energy 
than the AFM state but now only by 4(2 — y/3)t 2 /A per cluster. This is basically the reason 
for the significant reduction in the stability energy of the FM state against the AFM state 
in the isolated cluster as shown in Table I. 

Now we give some additional comments on the role of edge Bz in the stabilization of 
AFM state in isolated cluster. If there is no electronic coupling among the Bz's, the energy 
level shift of the edge pS state will be simply — t' 2 /A. Then the stability energy of the FM 
state against the AFM state for the isolated cluster will become At 2 / A, just half of that for 
the cyclic boundary condition, as is naively expected. However, because the pS states of 
the three Bz's form a coupled degenerate system, all the three pS states contribute to the 
polarization of the edge Bz and the actual energy shift is much larger as shown in Fig.5. 
This can be seen in the AFM DOS of Fig. 6. The state around -1.0 eV is the magnetically 
polarized Bz p5 state which consists mainly of the edge Bz state and also partly of the 
middle Bz state. Although the FM state is still more stable than the AFM state for any 
values of t' and t, the energy difference is reduced to 2{t/t') 2 (t 2 / A) when \t/t'\ « 1.0. 

In order to make the model analysis more quantitative, the parameters in the model are 
adjusted to reproduce the energy diagram obtained by the first-principles calculation. The 
total and partial DOS for the AFM and FM V 2 Bz 3 calculated with [7=3.0 eV are plotted 
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FIG. 5: (Color online) Energy diagrams for AFM and FM V2BZ3 cluster and (VBz)2 with cyclic 
boundary condition. e p = = 6q is assumed and the energy levels are given only in the lowest 
order of t 2 /A and t' 2 /A. 



TABLE II: Magnetic configuration (Arrow indicates the direction of local moment on V atoms and 
vertical lines represent Bz's) of V3BZ4, number of Bz's sandwiched with antiparallel V magnetic 
moments (represented by vertical dotted lines) and the energy difference (in eV) relative to the 
FM state from model calculation. Those in parentheses are from first-principles GGA+U calcula- 
tions with [7=3.0 eV. The structure is optimized for each magnetic configuration in the GGA+?7 
calculations, while t = t' =1.648 eV is used in the model analysis. 



magnetic configuration 


No. of | : i pairs 


Energy difference relative to FM state 


1 T 1 T 1 T 1 
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0.033 (0.018) 
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0.060 (0.037) 
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FIG. 6: (Color online) Total and partial DOS for AFM and FM V2BZ3 cluster calculated with 
GGA+/7 and [7=3.0 eV. 



TABLE III: Magnetic configurations (Arrow indicates the direction of local moment on V atoms 
and vertical lines represent Bz's) of V4BZ5, number of Bz's sandwiched with antiparallel V mag- 
netic moments (represented by vertical dotted lines) and the energy difference (in eV) relative 
to the FM state from model calculation. Those in parentheses are from first-principles GGA+C7 
calculations with U =3.0 eV. The structure is optimized for each magnetic configuration in the 
GGA+U calculations, while t = t' =1.648 eV is used in the model analysis. 



magnetic configuration 


No. of f : I pairs 


Energy difference relative to FM state 


1 T 1 T 1 T 1 T 1 





0.0 


ItltltUl 


1 


0.037 (0.036) 


itituui 


1 


0.051 (0.056) 


1 T : i U i T 1 


2 


0.073 (0.065) 


1 T 1 T : X : T 1 


2 


0.077 (0.082) 


1 T : i : T : i 1 


3 


0.105 (0.103) 
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FIG. 7: (Color online) Site % and cluster size n dependence of the energy cost AE(i,n) to have an 
antiparallel magnetic moment pair. The site i corresponds to the location of the Bz sandwiched 
with antiparallel magnetic moments of V and the cluster size n corresponds to V n Bz n+ i. t = t' is 
assumed. In the inset, the configuration is specified with i=l and n=5. 

in Fig. 6. Without losing generality, e p can be taken as 0.0 to be the energy reference 
and t' is first assumed to be equal to t. t, e& and A are adjusted using the occupied four 
levels in the FM state and the values for the parameters are t=1.648 eV, = —0.650 eV 
and A=0.916 eV, respectively. Then in order to clarify the effects of variation of t'/t in the 
relative stability between FM and AFM states, we assume the following distance dependence 
of the p — d hybridization, t,t' oc d~ 4 with d the distance between C and V and readjust 
t and if keeping e^, A and the total energy of the FM state unchanged. According to our 
structural optimization for [7=3.0 eV, the V-C distances are 2.304 A for the middle Bz and 
2.248 A for the edge Bz. The final values of t and t' are 1.555 eV and 1.716 eV, respectively. 
By using these parameters, we find that FM V2BZ3 is lower than the AFM one by 0.018 
eV. The corresponding value for the case of t — i' is 0.023 eV. For the cyclic case, FM 
state is more stable by 0.245 eV. These results are qualitatively consistent with our first- 
principles calculations. If a small contribution from da states to AFM coupling may be 
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taken into account, agreement between the model and the first-principles calculation may 
become better. 

The model analysis can be easily extended to larger clusters to study the size effect. Here 
we show some results for n=3 and 4 cases. In larger clusters there are several magnetic 
configurations. Table II and Table III summarize the energy cost needed to reverse the spin 
magnetic moment of V from the ground state for n=3 and 4, respectively. In the left-end 
column of the tables, arrows indicate the direction of V magnetic moment, while the vertical 
line corresponds to Bz. Those Bz's sandwiched with antiparallel V magnetic moments are 
denoted by dotted vertical lines. First, we note that in the GGA+U scheme, the ground 
state is FM for both n =3 and 4 being consistent with other GGA-BLYP 18 calculations, 5 
though one of the GGA-PBE calculations claims that the AFM state is the ground state for 
n =3.- Second, the energy cost for reversing the magnetic moment of one of the edge V's 
increases with n. This trend is seen in both the GGA+U calculation and the model analysis. 
Third, the energy cost is basically controlled by the number of dotted vertical lines , i.e., 
the number of antiparallel pairs, in the magnetic configuration. Fourth, between the two 
configurations with one dotted vertical line (the 2nd and 3rd configurations in Table III), 
the one with antiparallel spin configuration at the cluster edge costs less energy. The second 
and the fourth aspects are clearly seen in Fig. 7 where the results by model calculations are 
also shown for larger clusters. 



IV. CONCLUSION 



By examining the electronic structures of the infinite chain of VBz obtained with the 
GGA+U calculations, we first found that the orbitals with different symmetries play distinct 
roles in the physical properties of the VBz system. The da orbital of V, which has very weak 
hybridization with pa orbital of Bz, is localized and singly occupied with a large exchange 
splitting. It contributes one \1b magnetic moment to the whole system but plays minor role in 
the magnetic coupling between neighboring V local moments. Perhaps, it may contribute to 
weak anti-ferromagnetic superexchange. The dn and pn states hybridize strongly to push the 
7T molecular orbitals far away from the Fermi level either fully occupied or empty. Therefore, 
although they contribute to the stability of the multiple-decker sandwich structure, they are 
inactive in magnetic or electronic properties. The hybridization between dS and pS states 
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is medium among the three types of bonds and these states are crucially important in the 
magnetic coupling. We pointed out that the mechanism of ferromagnetism stability proposed 
by Kanamori and Terakura- 3 is directly applicable to the ferromagnetism of the infinite 
chain of VBz. We have also found that the presence of edge Bz in V2BZ3 cluster introduces 
an interesting modification to the mechanism. The edge Bz's are not simply inactive in the 
magnetic coupling but tend to favor the antiferromagnetic coupling of V magnetic moments. 
We set up a simple tight-binding model to analyze the problem in detail. Finally we studied 
the size dependence of magnetic interaction in the V n Bz n+ i clusters. 
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